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Abstract: Recent data on /e~ and p cosmic rays suggest that dark matter annihilate 
into the standard model (SM) particles through new leptophilic interaction. In this paper, 
we consider a standard model extension with the gauged U{1)L^-Lr group, with a new 
Dirac fermion charged under this U{1) as a dark matter. We study the muon {g — 2)^, 
thermal relic density of the cold dark matter, and the collider signatures of this model. Z 
productions at the Tevatron or the LHC could be easily order of 0(1) — 0(10'^) fb. 
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1. Introduction 

The Standard Model (SM) should be extended in order to accommodate the nonbaryonic 
cold dark matter (CDM) of the universe. There are many options for cold dark matter 
candidates in particle physics models: axion, neutralino or gravitino lightest supersymmet- 
ric particle (LSP) fl], ^, lightest Kaluza-Klein particle (LKP) and lightest particle in 
a hidden sector Q, to name a few. However, we do not have enough information on the 
detailed nature of CDM, except that ^cDAih'^ = 0.106 ± 0.008 More information from 
direct and indirect dark matter searches and colliders are indispensable for us to diagnose 
the particle identity of the CDM, namely its mass and spin and other internal quantum 
number (s). 

Recently, PAMELA reported a sharp increase of positron fraction e"*"/ (e"*" + e~) in the 
cosmic radiation for the energy range 10 CeV to 100 GeV and no excess in p/p Q 
from the theoretical calculations. And also very recently, Fermi-LAT ||7| and HESS Q 
data showed clear excess of the e"*" + e~ spectra in the multi-hundred GeV range above the 
conventional model although they do not confirm the previous ATIC [10| peak. 



These data (at least a part of them) could be due to some astrophysical origins such 
as pulsars H, 0]. A more exciting possibility from particle physics point of view would 
be interpreting them as indirect signatures of cold dark matters through their pair an- 
nihilations or decays. In this paper we take the second avenue, namely particle physics 
explanations of positron excess. Model independent study|13| show that multi-TeV scale 



- 1 - 



DMs dominantly annihilating into SM leptons, especially into r'^r" or 4 /i's, are most 
favored. If one assumes the standard cosmology, large boost factor (BF) of O(IO^) is also 
required to enhance the event rates. 

From the view point of model building and grand unification, it is non-trivial to con- 
struct this kind of leptophilic model. However it is still phenomenologically viable, and 
we have to verify or falsify this class of models by comparing the predictions with the 
data. In this paper we explicitly construct a leptophilic model and work out the physical 
consequences in detail. There are already many papers available studying the implications 



of the PAMELA/FERMI data in different models and/or context [|14|. 



The simplest model for the leptophilic (or hadrophobic) gauge interaction is to gauge 



the global U{1)l Lt symmetry of the standard model (SM), which is anomaly free |15 



[T6| , |T8|. Within the SM, there are four global U{1) symmetries which are anomaly free: 

Le ~ Lfj^, — Lt, Lt — Lg, B — L 

One of these can be implemented to a local symmetry without anomaly. The most popular 
is the U{\)b-l, which can be easily implemented to grand unified theory. Two other sym- 
metry involving Lg are tightly constrained by low energy and collider data. On the other 
hand, the — symmetry is not so tightly constrained, and detailed phenomenological 
study has not been available yet. Only the muon [g — 2)^ and the phenomenology at muon 
colliders have been discussed This model can be extended by introducing three 

right-handed neutrinos and generate the neutrino masses and mixings via seesaw mecha- 
nism [jl^. Also U{1)L^-Lr can be embedded into a horizontal SU{2)h acting on three 
lepton generations, which may be related with some grand unification. 

In this paper, we extend the existing U{1)l^-l^ model by including a complex scalar 
(j) and a spin- 1/2 Dirac fermion tp£,, with U{1)l^-l^ charge 1. There is no anomaly 
regenerated in this case, since we introduced a vectorlike fermion. The complex scalar (j) 
gives a mass to the extra Z' by ordinary Higgs mechanism. And the Dirac fermion plays 
a role of the cold dark matter, whose pair annihilation into fi or t explains the excess of e"*" 
and no p excess as reported by PAMELA |5|, |6| . Then we study CDM cosmology (thermal 
relic density and (in)direct signatures) and collider phenomenology of the C/(1)l^_l^ model 
with Dirac fermion dark matter in detail. 

In Sec. 2, we define the model and discuss the muon — in our model. In Sec. 3, we 
calculate the thermal relic density of the CDM ipu, and identify the parameter region that is 
consistent with the data from cosmological observations. We also present the signatures for 
indirect search experiments: e~^e~ , neutrinos, and gamma rays from the DM annihilations, 
including the Sommerfeld enhancement. In Sec. 4, we study the collider signatures of 
the model at various colliders encompassing Tevatron, B factories, LEP(2), the Z° pole 
and LHC, including productions and decays of Z' , the SM Higgs boson and the newly 
introduced U{1)l^-l^ charged scalar boson. Our results are summarized in Sec. 5. We 
note that this model was discussed briefly by Cirelli, et.al in Ref. |]l^] in the context of the 
muon {g — 2)^ and the relic density. In this paper, we present the quantitative analysis on 
these subjects in detail, as well as study various signatures at various colliders including 
the Tevatron and the LHC. 
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2. Model and the muon {g — 2)^ 

The new gauge symmetry U{1)l^-l^ affects only the 2nd and the 3rd generations of leptons. 
We assume generation index) carry Y' = 1(— 1). We further introduce 

a complex scalar (p with (1, 1, 0)(1) and a Dirac fermion '(/'£) with (1, 1, 0)(1), where the first 
and the second parentheses represent the SM and the U{1)L^-Lr quantum numbers of 4> 
and tpD, respectively. The covariant derivative is defined as 

D^ = d^ + ieQA^ + i^—ih - slrQ)Z^ + igY z' (2.1) 
swcw 

The model lagrangian is given by ^ 

■^Model = -^SM + -^New (2.2) 

+D,,cP*Df'cP - iil<p*4> - A<^(0»2 - \H^(t>*(t)H^H. (2.3) 

In general, we have to include renormalizable kinetic mixing term for U{\)y and U{\)l^^-Lt 
gauge fields, B^i,Z''^'^, which will lead to the mixing between Z and Z' . Then the dark 
matter pair can annihilate into quarks through Z — Z' mixing in our case, and the p flux will 
be somewhat enhanced, depending on the size the Z — Z mixing. However, electroweak 
precision data and collider experiments give a strong constraint on the possible mixing 
parameter, since the mixing induces the Z' coupling to the quark sector. Furthermore, if 
one assumes that the new U{1)l^^l^ is embedded into a nonabelian gauge group such as 
SU{2)h or SU{3)h, then the kinetic mixing term is forbidden by this nonabelian gauge 
symmetry [p!6| . In this paper, we will assume that the kinetic mixing is zero to simplify the 
discussion and to maximize the contrast between the positron and the antiproton fluxes 
from the dark matter annihilations. ^ 

In this model, there are two phases for the extra U{1)L^-Lr gauge symmetry depending 
on the sign of /i^ : 

• Unbroken phase: exact with ((/>) = 0, /x^ > and M^i = 0, 

• Spontaneously broken phase: by ^u^ < 0, nonzero (i;^) = ^ 0, and M^/ / 

In the unbroken phase, the massless Z' contribute to the muon (g — 2)^ as in QED up to 
the overall coupling: 

Aa^ = |-. (2.4) 

Currently there is about 3.4(T difference between the BNL data [^2| and the SM predic- 
tions |23| in {g — 2)^: 

Aa^ = a;fP - af^ = (302 ± 88) x IQ-^^ (2.5) 

^Similar idea for the DM was considered in po| pl[ in the context of Stueckelberg U{l)x extension of 
the SM model. 

^Because of this simplification, the direct detection rate from the CDM (in)elastic scattering off nuclei 
vanishes identically. However this is no longer true if the kinetic mixing is included. This case is discussed 
in brief in Sec. 3.2. 
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The Aa^ in (|2.4| ) can explain this discrepancy, if a' ~ 2 x 10^^. However, this couphng 
is too small for the thermal relic density to satisfy the WMAP data. The resulting relic 
density is too high by a several orders of magnitude. Also the collider signatures will be 
highly suppressed. Therefore we do not consider this possibility any further, and consider 
the massive Z' case (broken phase) in the following. 

In the broken phase, it is straightforward to calculate the Z' contribution to Aa^. We 
use the result obtained in Ref. [18|: 



a 2mf^x^{l-x) a 2^2 

^"'^ = 2^ io '^'^x2m2 + (l-x)M2, ^ 2^3Mf ^^'^^ 



The second approximate formula holds for <C M^'. In Fig. 1, shown in the blue band 
is the allowed region of M^/ and a which is consistent with the BNL data on the muon 
{g — 2)^ within 3 a range. There is an ample parameter space where the discrepancy 
between the BNL data and the SM prediction can be explained within the model. 



3. Dark matter : Relic density and (In) direct signatures 
3.1 Thermal relic density 

In our model, the Dirac fermion ipjy and its antiparticle 'i/'D CDM candidates. The 
thermal relic density oi ipo and tp^ is achieved through the DM annihilations into muon, 
tau leptons or their neutrinos through s-channel Z'-exchange. They can also annihilate 
into the real Z' pairs when kinematically allowed. 

iJoi^D —> z'* I'^r ^vivi {I = fj.,T), 

i^DtpD ^ Z'Z'. (3.1) 



We modified the micrOMEGAs [24| in order to calculate the relic density of the U{1)l^-l^ 
charged V'D CDM. It is easy to fulfil the WMAP data on f^cDM for a wide range of the DM 
mass, as shown in Fig. [l|. The black curves represent constant contours of O/i^ = 0.106 
in the {Mz' , a) -plane for M^^ = 10, 100, 1000 GeV (from below). We can clearly see the 
s— channel resonance effect of Z' tpD'ipD near Mz' ~ 2M^^. The blue band is the 
allowed region by the {g — 2)^ at the 3 a level. We also show the contours for the Z' 
production cross sections at various colliders: B factories (Ifb, red dotted), Tevatron (lOfb, 
green dot-dashed), LEP(10fb, pink dotted), LEP2(10fb, orange dotted) and LHC (1 fb, 
10 fb and 100 fb in blue dashed curves). The cross sections in the parentheses except the 
LHC case roughly correspond to the upper bounds that each machine gives. Therefore the 
left-hand sides of each curve is ruled out by the current collider data. Note that a larger 
parameter space can be accessed by the LHC. These issues and other collider siugnatures 
are covered in the next section. 

The current experimental mass bound of SM-like Z' is 923 GeV from the search for 
a narrow resonance in electron-positron events [p5||. We emphasize, however, that in our 



model the Z' boson as light as ~ 10 GeV is still allowed by present data from various 
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colliders. It is mainly because the production cross section at the Tevatron is suppressed 
since Z' should be produced from the couplings to the 2nd and 3rd family leptons. 

In the range 100 GeV < M^^ < 10 TeV, a > lO'^ and 100 GeV < Mz' < 1 TeV, the 
relic density and Aa^ constraints can be easily satisfied simultaneously while escaping the 
current collider searches. We note that if the {g — 2)^ constraint is not considered seriously 
or if we assume there are other sector which saturate the (^f — 2)^ upper bound, then all 
the region in the right-hand side of the blue band is also allowed. 




logjoCM^ /GeV) 

Figure 1: The relic density of CDM (black), the muon (g — 2)^ (blue band), the production cross 
section at B factories (1 fb, red dotted), Tevatron (10 fb, green dotdashed), LEP (10 fb, pink 
dotted), LEP2 (10 fb, orange dotted), LHC (1 fb, 10 fb, 100 fb, blue dashed) and the Z° decay 
width (2.5 xlO~^ GeV, brown dotted) in the (logj^g ilogio-^z') plane. For the relic density, we 
show three contours with nh'^ = 0.106 for A/^„ = 10 GeV, 100 GeV and 1000 GeV. The blue band 
is allowed by Aa^ = (302 ± 88) x 10"" within 3 cr. 



3.2 Direct detection rates 

Since we ignored the kinetic mixing between the new U{1) gauge boson and the SM U{1)y 
gauge boson B^, there would be no signal in direct DM detection experiments in this 
model. The messenger Z' does not interact with electron, quarks or gluons inside nucleus. 
Also there would be no excess in the antiproton flux in cosmic rays in this case, while one 
could have an excess in the positron signal in a manner consistent with the PAMEL A/Fermi 
data. However there would be a small kinetic mixing between two U{1) gauge field strength 
tensor. If we assume a small kinetic mixing 10~^ = 10~^) between the Z' and photon, 
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then spin-independent cross section for direct detection rate will be given by 



where Z and A are the atomic number and the mass number of a nucleus. The couplings 
Ap(„) 's of the CDM to proton and neutron are given by 



eg 



" 2M2, 



^^^ + 4smZ:os.., (1-^-^^-) 



(3.3) 



" " ^8M2,sin^^cos0^ ^ ^ ^ 

where the upper (lower) sign corresponds to ipDii^n) scattering. Note that Xp and 
Xn are dominated by the — 7 and Z' — Z mixing, respectively. This is because photon 
couples to the nucleon charge, whereas the Z^ couples to the neutral current weak charge 
of a nucleon. The Z^ coupling to a proton is proportional to (1 — Asin^Ow), and thus 
highly suppressed compared to the coupling to a neutron. Also let us note that the 
cross section on a nucleus could be small, if there is a cancellation between Xp and Xn terms 
depending on the sign of 9^'^ and O^'z ^(i-s (3.3)-(3.5). 

If we consider a dark matter scattering on single proton target, one has 



4- , . ^ - . :„-«c™' (l^))' (^\ (^)^ ,3,5, 



which is dominated by Z — 7 mixing. The resulting cross section is close to the current 



upper bounds from XENONIO and CDMS [27| experiments. We have assumed that 
the DM V'D is much heavier than proton. Similarly, the SI cross section on a neutron target 
is given by 

"IL^ . 4.6 X io-«e„.^ /imy\' 

The scattering cross section on the proton hy Z — Z' mixing is suppressed by (1 — 4 sin^ ^u;)^ 
relative to the scattering cross section on the neutron target, and thus negligible. In either 
case, one can evade the bounds from XENONIO and CDMSII by taking a heavier Z' mass, 
smaller coupling a or smaller mixing ang le 61. If Z is Ught, one may have too large SI 
cross section, in conflict with XENONIO and CDMSII. 




3.3 Sommerfeld enhancement and boost factor (BF) 

The DM annihilation cross section at the freeze-out temperature is typically (av) ~ 3 x 
10~^^ cm^/sec. To explain the PAMELA or Fermi data, (av) ~ 3 x 10~^^ cm^/sec is 
required. Therefore we need a boost factor (BF) of order of 10^. A large BF can come from 
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the so-called Sommerfeld enhancement in the 5— wave DM annihilation. The Sommerfeld 
factor, given by the ratio of the radial wavefunction at infinity to that at the origin, 



Sk 



Xfc(oo) 



(3.7) 



can be calculated by solving the radial S— wave Schrodinger equation |28| with the attrac- 
tive Yukawa potential in our model: 



1 



Xkir) 



fc2 



-Xkir)- 



(3.8) 



Here k = M^^v and v is the relative velocity of two annihilating DM particles. The 
boundary condition for the above Schrodinger equation is 



Xfc(r) ^ ifcXA:('^) as r ^ oo, Xfc(O) = 1. 



(3.9) 



Fig. ^ shows the prediction for the Sommerfeld enhancement factor in our model for various 
values of DM: M^^ = 10, 100, 1000, 2000 GeV. For a given values of the DM mass and the 
M^/, the a' are chosen in such a way that they satisfy the relic density, i.e. each lines are 
predictions of the Sommerfeld enhancement factor along the constant relic density contours 
in Fig. 2. We can see that it is easy to get the enhancement factor ~ 10^. 




12 3 4 

logio(M2/GeV) 



Figure 2: The predictions of Sommerfeld enhancement factor in our model along the constant 
contours in Fig. 2. 
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3.4 Indirect signatures: positron, neutrino and photon fluxes 

Now we show the prediction of jie^ + e^) spectra in our model in Fig. ^ The positron 
flux at the Earth can be calculated from the solution of the diffusion equation as |2£] 

^e+(0,e) = desfies)iiXD), (3.10) 

where represent the position of the sun, e = Ef,+ /EQ {Eq = 1 GeV). The function f{es) 
is the positron energy spectrum from the dark matter annihilation at the source. The 
"diffusive halo function" I(Ad) with A|, = 4KETEi€^~^ - e|"^)/(l - 5) encodes the in- 
formation on the propagation of positron from the source to the Earth. The parameter 
K = r]{av) (pq/M^j^)'^ is a factor relevant to the particle physics [77 = 1/2 (1/4) for Ma- 
jorana (Dirac) DM]. The other parameters are: the speed of the positron /3g+ and and 
TE = 10^^ sec. We used the NEW DM density profile @: 

(/3-7)/" 

(3.11) 



frey ^ l + (re/r,)" V^-^)/" 



where (a,/3, 7,rs) = (1,2,1,20 kpc) and Pq = 0.3 GeV/cm'^ is the DM density near the 



Sun. To obtained the halo funtion /(A^)), we used the method suggested in |31]. We also 
used the cosmic ray propagation parameters which correspond to the medium primary 
antiproton fluxes |2|: 5 = 0.70, Kq = 0.0112kpc2/Myr, and L = 4kpc. Eor the plot, 
Eig. ^, we fixed the DM mass M^^ = 2 TeV. The required BE is about 5200, which is a 
little bit larger than the maximal Sommerfeld enhancement in Eig. ^ can give. However, 
an additional enhancement factor of about 2-3 can be easily obtained from the clumpy 
structure of dark matter density. We also have checked with other DM masses ranging 
from 1 TeV to 3 TeV, which can also fit the PAMELA data very well. 



Although the PAMELA data alone can be fitted with a wide range of the DM mass [13|, 
the simultaneous fit including the Eermi and HESS data is non-trivial and gives more strong 
constraint on the DM mass. Eig. ^ shows a fit to the PAMELA, Eermi and HESS data in 
our model. Eirst, we obtained the absolute positron fiux from the PAMELA data and the 
known background positron and electron spectrum Q. The resulting spectra agree with 
the Eermi data when we rescale the background by a factor r ~ 0.7. This rescaling also 
makes all the Eermi data lie above the background. The HESS data has also been allowed 
a rescale factor rn ■ Now we fitted the three parameters r, and the BE to the data, 
assuming the data are independent with each other. Eor the DM mass M^^ = 2 TeV, we 
obtained an excellent fit xl^iJd-o.f = 53/50, r = 0.7, = 0.84 and BE = 5200. Eor the 
lighter and heavier DM masses, the fit quality becomes worse. We have also checked that 
the isothermal DM profile also gives similar results. The reason can be traced back to the 
fact that the positrons we observe comes mainly from the sources not that far from the Sun 
where the DM density does not differ much for different DM profiles, although the NEW 
profile is far more cuspy than the isothermal profile when approaching the Galactic center. 

Since the DM pair annihilates into the 2nd and 3rd generation leptons including neu- 
trinos in our model, we expect large neutrino fiux. The neutrino flux can be detected 
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I PAMELA (2009) 

Background 

^ DM M=2000 (GeV) 




Figure 3: The fit to the PAMELA data in our model. 



^ DM M=2000 (GeV) 
• Fermi (2009) 
■ PAMELA (2009) 
HESS(09) X 0.84 




Figure 4: The fit to the PAMELA, Fermi and HESS data in our model. 

through the upward-going muons in the Super-Kamiokande (SK). Also the neutral pion 
from the tau decay can produce sizable photon flux, which can also be compared with the 
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existing gamma-ray searches such as HESS^. The main contribution to the neutrino and 
the gamma-ray flux comes from the Galactic center where the DM density is the highest. 
For this reason we consider only neutrinos and gamma-rays from the Galactic center. The 
differential fluxes of neutrinos and photons from the Galactic center can be easily calculated 
from the formula |^] for the case of Dirac DM: 

= M^TTT^T^ / dsp\r{s,^P)), (3.12) 

dEi dE levrM^^ Jx^r^e of sight 

where z = z/, 7 and s is the distance from the Earth in the angular direction ip from the 
line connecting the Earth and the Galactic center (GC). 

The neutrinos from the GC can be detected at the superkamiokande (SK) as the muon 
neutrios transform into the muons through the weak interactions in the rocks below the 
SK. The neutrino-induced muon flux is written as 

r dF 

Vm- = j ^^-M^/(^-.)- (3-13) 

The function f{E^^) is the probability of a muon neutrino with energy E^^ transforming 
into muon with energy larger than Eth, and is given by [^] 

where da,y^(p^^p(^n)~>^i(ji)x / dE^^ is the scattering cross section of a neutrino with proton 



(neutron) to create a muon with energy Ef^ |35]: 



d^. + 

dE, vr 1^ ElJ ^ > 

For the number density of proton (neutron) in the rock, we use rip"'^^^ = nn°^^^ = 
2.65A^A/2cm-3 {Na = 6.022 x 10^3). R{E,,Eth) is the distance a muon with E, can 



travel inside the rock before losing energy below Eth, and is fitted to be [34 1 

R{E,, Eth = lOGeV) = W+^y+^^y^ (km), (3.16) 

where y = logio(£^/./lGeV), a = -3.29186,6 = 1.52594, and c = -0.147224. 

Fig. ^ shows the predictions for the neutrino-induced muon flux for the DM masses 
M^^ = 3,2, 1.5, 1 TeV (from above). We obtained the annihilation cross section in such 
a way that each DM mass fits the PAMELA, Fermi and HESS data as described above. 
We used the NEW (solid red curves) and the isothermal (dashed blue curves) profiles for 
the plot. We can see that the 3 TeV DM is already ruled out by the SK bound because it 
needs too large BF. The 2 TeV DM which fits the CR data best is only marginally allowed. 



^See also |33] for the gamma-ray constraint. 
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Cone half angle from GC [deg] 

Figure 5: Thick solid red curves (thick dashed blue curves) are predictions of the neutrino- induced 
up-going niuon flux from the annihilation of dark matter with masses 3, 2, 1.5, 1 TeV from above, 
for the NFW (isothermal) dark matter profile. The thin solid line is the superkamiokande bound. 



The lower DMs are allowed with the NFW profile. However, if the isothermal profile is 
used, all the DM are allowed because this profile is flat near the Galactic center and the 
neutrinos are not much produced. 

Fig. P shows the predictions for the gamma-ray flux from the Galactic center (0.1° 



region from the GC) |36|] and the Galactic Center ridge < 0.3°, |/I < 0.8°) |37|1. We can 
see that the constraints on the DM annihilation for the NFW profile become more severe 
than in the neutrino case. That is the NFW predicts too much gamma-ray, exceeding 
even the current data for the massive DM. However, if more flat profile like the isothermal 
profile is used, the predictions are below the current data. 



4. Collider Signatures 

New particles in this model are Z , s (the modulus of 0) and V'D- Z couples only to muon, 
tau or their neutrinos, or the C/(1)L^-Lr charged dark matter. The new scalar s can mix 
with the SM Higgs boson /isMj affecting the standard Higgs phenomenology. 

Let us discuss first the decay of Z' gauge boson and its productions at various colliders. 
In the broken phase with M^/ ^ 0, Z' can decay through the following channels: 

Z —)■ fi'^ fi~ , t'^ T~ , Uai^a (with a = fi OT r), ^jdiPd j 

if they are kinematically allowed. Since these decays occur through U{1)l^-l^ gauge 
interaction, the branching ratios are completely fixed once particle masses are specified. In 
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Galactic Center 



Galactic Ridge |b|<0.3°, |1|<0.8° 




E [GeV] E [GeV] 



Figure 6: The gamma ray flux from the GC (left panel) and GC ridge (right panel). Thick solid 
red curves (thick dashed blue curves) are predictions of the gamma ray flux from the annihilation 
of dark matter with masses 3, 2, 1.5, 1 TeV from above, for the NFW (isothermal) dark matter 
profile. 



particular, 

T{Z' ^ fi+fi-) = r(z' ^ r+T-) = 2r(z' ^ Uf^Uf,) = 2T{Z' ^ u^Ur) = T{z' ^ VdV^z?) 
if M^' ^ m^, niT, M^m- The total decay rate of Z' is approximately given by 

TMZ') = ^ M^, X 4(3) « 1^ GeV ( ( 

if the channel Z' — > ipn'ipD is open (or closed). Therefore Z' will decay immediately inside 
the detector for a reasonable range of a and M^/ . 

Z' can be produced at a muon collider as resonances in the /^/x or rr channel |18] via 

The LHC can also observe the Z' which gives the right amount of the relic density as can 
be seen in Fig. 1. Its signal is the excess of multi-muon (tau) events without the excess of 
multi-e events. 

The dominant mechanisms of Z productions at available colliders are 

qq (or e'^e") — > 7*,Z* — > /i^/i^Z ,t^t^Z 

^ Z* ^ U^U^Z , VrVrZ 

There are also vector boson fusion processes such as 

W'^W- u^Vf,z' (or /i+^"Z'), etc. 



Z^ Z^ — > v^v^^Z (or [i^ [I Z ), etc. 



W^Z'' v^liZ (or /i+^"Z ), etc. 

and the channels with fi ^ t. We will ignore the vector boson fusion channels in this paper, 
since their contributions are expected to be subdominant to the qq or e~^e~ annihilations. 
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In Fig. 1, we present the Z' production cross sections at B factories, pole, LEP(2), 
Tevatron and LHC. We find that the hght M^i region that can accommodate the muon 
((7 — 2)^ is almost excluded by the current collider data. The remaining region can be 
covered at the LHC with high integrated luminosity > 50 fb^^. 

The signatures of Z' will be an s— channel resonance in the dimuon invariant mass 
spectrum, or its deviation from the SM predictions as in Drell-Yan production of the muon 
pair. Therefore one could expect that the number of multi-muon events at colliders is 
enhanced compared with the SM predictions. The e^e~ channel will be diluted compared 
with the z^"*"/^" channel, since the e^e~ final state in the Z decay can appear only through 
the Z' decay into a tau pair and r evv in this model. 

Now let us discuss the Higgs phenomenology in our model. In general, there can be 
a mixing between the SM Higgs boson /igM and a new scalar s (the modulus of 0) due to 
the coupling in ( |2.3D . As a consequence, the Higgs searches at colliders can be quite 
different from those of the SM. For example, one can imagine 

99 ^SM s* ^ Z Z\ 

followed by Z ^ fj.'^ fi~ , t~ , v^v^, Vrh'r,'>pD'>pD- This makes an additional contribution to 
the Z' production at the LHC. However, we did not include this Z' pair production through 
gluon fusion in Fig. 1 for simplicity, since it depends on the unknown free parameter \H(f>, 
and thus is more model dependent. In any case, the generic collider signatures of the 
new Z' are the excess of multi-muon or tau events, compared with the SM. It is strongly 
desirable to search for /U/x/i/i, tttt, or fifirr or large missing Et from fifivu events at LEP, 
LEP2, Tevatron and at the LHC. 

We can introduce an angle f3 so that the ratio of two scalar VEV's is given by tan /3 = 
V(f)/vhg-^, and an angle a that parametrizes the mixing of /igM and s: 

hsM = Hi cos a — H2 sin a, 

s = i?i sin a + cos a, (4-1) 

where -^1(2) denotes the lighter (heavier) mass eigenstate of two scalars. 

In Fig. 0, we show the branching ratios (BRs) of the two-body decay modes of Hi^2 
for Mz' = 300 GeV. We have fixed Mh^ = 700 GeV {Mh^ = 150 GeV) for the plots of the 
Hi {H2) decay [ the left (right) column ]. Note that the modes Hi 2 Z' Z' (solid blue) 
and H2 — > HiHi (solid green), which are absent in the SM Higgs decay, can dominate for 
large a and small tan/3. If HiS and Z' are heavy enough compared with the CDM in our 
model, a decay Hi — > Z' Z' followed by one or both of the Z' decaying into a pair of CDM 
or a pair of neutrinos could occur. Therefore the Higgs could have somewhat large invisible 
branching ratio, compared with the SM Higgs boson. Therefore Higgs signatures at the 
Tevatron or the LHC could be quite exotic . 

5. Conclusions 

Recent possible anomalies in the cosmic ray data reported by PAMELA, Fermi-LAT and 
HESS Collaborations may be due to astrophysical origins such as pulsars. However, it is 
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Figure 7: In the left (right) column are shown the branching ratios of the lighter (heavier) Higgs 
i?i(2) into two particles in the final states: ti (solid in red), bb (dashed red), cc (dotted red), ss 
(dot-dashed red), rf (solid orange), fifl (dashed orange), WW (dashed blue), ZZ (dotted blue) and 
Z Z (solid blue) for difference values of the mixing angle a and tan/3. We fixed Mz' = 300 GeV. 
We also fixed Mh^ = 700 GeV {Mh, = 150 GeV) for the plots of the left (right) column. 



also tantalizing to consider them as the first hint for the existence of weakly interacting cold 
dark matter. In this paper, we considered a leptophilic CDM model with extra ?7(1)l^-l^ 
gauge symmetry which is one of the anomaly free global symmetry in the SM. We have 
introduced a new complex scalar (p and Dirac fermion which are charged under the new 
U{1)l^-l^. The U{1)l^-l^ charged Dirac fermion ipjj can be a good CDM that might 
explain the positron excess reported by HEAT, PAMELA and FERMI, without producing 
excess in antiproton flux as observed by PAMELA. This model is constrained by the muon 
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(g — 2)^ and collider searches for a vector boson decaying into ^Ji^ ^~ at the Tevatron, 
LEP(2) and B factories. The collider constraints favors Vdm heavier than ~ 100 GeV. 
We calculated the relic density of the CDM with these constraints, and still find that the 
thermal relic density could be easily within the WMAP range. We also considered the 
production cross section of the new gauge boson Z' at the LHC, which could be 1 fb -1000 
fb. The new gauge boson Z' will decay into fijl, rf, their neutrino partners or even to a 
pair of CDM's. Therefore the final states will be rich in muons or taus, or missing Et- 
Most parameter space of this model is within the discovery range at the LHC with enough 
integrated luminosity > 50 fb~^. It is remained to be seen whether there are excess in the 
multimuon or multitau events at the Tevatron or at the LHC. 
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